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Abstract 

The superconducting charge-phase 'Quantronium' qubit is considered in order to develop a model 
for the measurement process used in the experiment of Vion et. al. [Science 296 886 (2002)]. 
For this model we propose a method for including the bias current in the read-out process in a 
fundamentally irreversible way, which to first order, is approximated by the Josephson junction 
tilted-washboard potential phenomenology. The decohering bias current is introduced in the form 
of a Lindblad operator and the Wigner function for the current biased read-out Josephson junction 
is derived and analyzed. During the read-out current pulse used in the Quantronium experiment 
we find that the coherence of the qubit initially prepared in a symmetric superposition state is lost 
at a time of 0.2 nanoseconds after the bias current pulse has been applied. A timescale which is 
much shorter than the experimental readout time. Additionally we look at the effect of Johnson- 
Nyquist noise with zero mean from the current source during the qubit manipulation and show 
that the decoherence due to the irreversible bias current description is an order of magnitude 
smaller than that found through adding noise to the reversible tilted washboard potential model. 
Our irreversible bias current model is also applicable to the persistent current based qubits where 
the state is measured according to its flux via a small inductance direct current superconducting 
quantum interference device (DC-SQUID). 

PACS numbers: 85.25.Cp, 74.50.+r, 03.65.Yz, 03.67.Lx, 
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I. INTRODUCTION 



Quantum computers and the quantum algorithms that run on them have been pro- 
posed as a technology to perform computational tasks not tractable with classical computer 
circuits^]. Recent experiments have provided significant advances towards developing the 
fundamental element of this technology, the quantum bit or qubit. So far, qubit systems 
based on nuclear magnetic resonance arid ion traps^,!^ have been used to show multiple 
qubit operation, whilst efficient linear optic quantum computing fij has been demonstrated 
with the successful operation of the two qubit controlled-not gate 7|. Experimental advances 
have also been made in solid state systems which utilise a wide variety of quantum effects 
in many different materials. The main attraction of solid state systems is the possibility to 
scale such technology using modern-day device fabrication techniques once the implemen- 
tation of component gates has been demonstrated. Promising solid state systems include 
the use of phosphor dopants in , optically 



controlled exciton systems [12] as well as a variety of systems based on the coherent electron 
state in superconducting materials ll 311 . 

nn 

iese superconducting systems the implementation of single qubit operation |14|, [15|, 
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181 ] . some with single shot readout, has been demonstrated. Also devices with a non- 
switchable inter-qubit interaction between two qubits have been shown 3 5 IsiJ , providing 
the initial evidence for a two-qubit entangled state in these structures. To ensure scalability 
to more complex configurations into the future there is a need to identify ways to develop 
more accurate gates, provide higher fidelity readout and ensure longer coherence times in 
the devices being developed. For instance the 'Quantronium' charge-phase qubit developed 
by Vion et a/.[l8| was designed to be insensitive to first order fluctuations in the external 
control parameters of the system provided that the control parameters for the device, in this 
case the voltage and applied flux, were used about an 'optimal point' of the system with 
this property. In this experiment the quality factor of quantum coherence Q for the device, 
defined as the number of elementary gate operations that could be performed before the 
device state decoheres, was found to be of the order 10 4 . 

In this paper we examine the readout process in the experiment of Vion et. al. through 
the Lindblad operator formalism 2lJ and we introduce the bias current into the model in 
a fundamentally irreversible way that acts to decohere the state of the qubit. Using this 
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method we implement a heuristic model for the measurement process that is induced by 
the application of the bias current to the quantronium circuit. This model allows for the 
bias current to 'count' the number of electrons that pass through the system during the 
measurement process and in doing so destroys the coherence between the different states 
of the system. Therefore to examine this model we are ignoring the typical terms which 
appear in the system master equation that describe the widely known forms of decoherence 
for the qubit through its coupling to the environment, such as the ohmic dissipation of the 
leads [22]. Our aim is to gain further insight into the role of the irreversible readout process 
and the decohering process associated with its operation. 

The irreversible dynamics arising from the current bias provides a decoherence mechanism 
that collapses the quantum superposition to a probabilistic mixture on a time scale shorter 
than the time for the state to tunnel out of the metastable qubit states into unbound states 
of the washboard potential and create a voltage on the read-out voltmeter. This means 
that the measurement of the system is performed during the application of the bias current 
before any classical information about the qubit state is returned to the experimentalist. 
Such 'measurement induced decoherence' is analogous to that discussed in semiconducting 
systems jQ]. In addition to this we also analyse the implications this irreversible current 
source has for the effect of Johnson-Nyquist noise from the current source during the qubit 
manipulation when the bias current has a zero mean and intended to be decoupled form the 
device. 



II. THE CURRENT BIASED JOSEPHSON JUNCTION 

The measurement process in superconducting qubit structures such as the Quantronium 
and the direct current superconducting quantum interference device (DC-SQUID) (which 
is used to measure the persistent current qubits|l3] and proposed to measure magnetic 
nanoparticlesj^j]) rely upon the transition of a Josephson junction based system from the 
superconducting state into the voltage state, where the information associated with the 
effective critical current of the device provides the quantum state measurement. The semi- 
classical model for a single Josephson junction is the one-dimensional analogy to a particle 
of mass (h/2e) 2 C moving along the 7 axis in the potential 2^1 

U( 1 )=E J (l-cos 1 )-E blj (1) 



where Ef, = Iuas^o/^- The term describes the slope of the washboard potential, 



which has been used widely in the quantum regime 



22 1 . For instance, it has been used 



to describe the escape rates of macroscopic tunnelling events in current biased Josephson 



junctions 



271 ] . The inclusion of the linear potential in Eq. (JTJ) to create the tilted 



washboard potential does not contribute any dephasing term to the dynamics and implies 
that the measurement process is intrinsically reversible. That is, by turning the current 
source on and then off again, the qubit is back in its initial state (provided a macroscopic 
quantum tunnelling event has not occurred). 

In this paper we propose an alternate description of the current bias in the Quantronium 
and other current biased systems such as the DC-SQUID, one which gives rise to the wash- 
board potential E^j term as well as intrinsically irreversible dynamics. This irreversiblity 
arises as a direct consequence of the measurement process, and the starting point for our 
model is the master equation 

P =4lH,p]+L P L<-^-^, (2) 

where for I^ as > we have defined L = a/ | E b | / hrf (see reference (28 1) and the charge- 
tunnelling non-unitary operator on the large Josephson junction is 77' \n) = \n + 1) and 
rj \n) = \n — 1). The state \n) represents the number of Cooper pairs that have tunnelled 
through the large Josephson junction (i.e an eigenstate of the Cooper pair number opera- 
tor N) and the cooper pair tunnelling operator 77 satisfies [N, 77^] = 77^ and [N, rj) = —rj. 
For Ibias < we have defined L = ^y\E b \/hr]. These Lindblad operators account for the 
movement of Cooper pairs across the Josephson junction at an average rate given by the 
current Iu as /2e. That is, the operators and 77 count the number of electrons added by 
the external bias current to the large Josephson junction at an average rate \Eb\/h. 

By introducing the Lindblad equation, given by Eq. (J2J), we are proposing a heuristic 
method to model the bias current which attempts to capture the notion that the current 
source counts the number of electrons tunnelling through the Josephson junction. In the 
remainder of this section we reconcile such a model by showing that it is in fact in agreement 
with a classical current biased Josephson junction and that the Lindblad terms contained 
in Eq. (J2J tend to destroy superpositions of different phase states. 
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Using Eq. (J2J), the master equation therefore reads 

. = ( -i[H, p] + Jfl (ijtpg - p) if J to > 

P |_*[ir,p] + Jfi(^t_p) if j Waa <o. 

In this equation the Hamiltonian H describes the Josephson junction or qubit dynamics 
but does not include the bias current washboard potential terms. Describing the current 
bias in superconducting circuits through this Lindblad superoperator is compatible with the 
phenomenology of the current biased Josephson junction in the classical limit (C — > oo) 
where the phase across the Josephson junction is fixed by the applied current. For instance, 
we can consider a single Josephson junction which is current biased (I bias > 0) and described 
by Eq. © where the Hamiltonian is given by 

ff=^V 2 + ^(l-co 87 ) 

and N is the Cooper pair number operator on the Josephson junction. In the steady state 
of this equation for the single Josephson junction, p = 0, we can compute the quantity 

~~jP~ ~ Tr(piV) — 0, 

to look at the role of the bias current in our model. Using the cyclic property of the trace 
we find that 

-i Tr(p [N, H}) + \E b \ Tr(r/ f p (N + 1) r] - pN) = 

and from the commutation relations for 7] and N, together with the definition of the current 
operator 

, 2tt dH T 

Ic sin 7 



$ <97 



we therefore show for Ibi as > that we have the expected result in the classical limit 
(C — > oo); that is 

(/) = Tr(p/) = Ibias. 

Also by considering an oppositely biased current (hias < 0) we find that by the inclusion of 
the Lindblad terms for the bias current in the master equation (Eq. (J3J)) we have (I) = Ibias 
and therefore retained the expected behaviour of the the bias current in the single Josephson 



junction system. That is, the current through the Josephson junction in our model is that 
applied by the current source. 

Additionally, in our model, the linear washboard term arises naturally from the 
Lindblad superoperator description of the bias current. By expanding the Lindblad super- 
operators in terms of its phase representation, then to first order, we can obtain the reversible 
dynamics of the washboard potential through the £^7 term. The higher-order terms from 
this expansion provide us with the intrinsic irreversible terms of the current source in our 
model. Hence having introduced the current source as an irreversible one, the system can 
be approximated by the washboard potential model of a current biased Josephson junction 
with an added irreversibility. For instance, by making an approximation to the full master 
equation (Eq. for the Quantronium circuit we can write the operators 77' and rj in their 
phase representation and approximate them to second order. That is we can write 

Tjt = e +i7 w 1 + ij - — (5) 

2 

and 

r ] = e- i ^\-i 1 - 1 -. (6) 

Under this approximation, and considering the cases for Iuas being positive and negative, 
the master equation of the system is 

* \n t? 1 1 \ Eb \ ( P^ 2 \ cr\ 

P = - l W-E blf ,p] + — ^ f pry- — -— y (7) 

Here we emphasise that the first order approximation to the operator L is the £^7 term which 
appears in the tilted washboard potential model. The additional three terms appearing at 
the end of the master equation are the irreversible decohering terms of this model under our 
second order approximation. 

By making the second order approximations (Eq. (JSJ) and Eq. (JBJ)) for the operators 
and 77, we have expanded them in terms of the operator 7 about the point (7) = 0; in this 
expansion we have used the small parameter A which is the variance of a sharply peaked 
Gaussian state in the phase representation. For instance if we consider only the decoherence 
term in the Iu as > master equation (Eq. Q) we have 

p=^(e^pe---p) = ^P[ 7 ]p. (8) 

The steady state of this master equation can be written as po = \4'o)(4 , o\ where j^o) is 
the sharply peaked Gaussian steady state wavefunction. This wavefunction results from 

7 



the small charging energy relative to the Josephson energy of the junction. Note that this 
wavef unction, tightly peaked around a given value of 7, is consistent with the Josephson 
relation for a classical current passing through a Josephson junction. Therefore we write the 
steady state wave function as 

where the variance A is small so that the wavefunction is sharply peaked in phase. Using 
this wavefunction to construct the steady state density matrix po we can approximate the 
term V[j}p in the master equation (Eq. (JHJ)) as follows: 

= - t ^A [-7, p ] - A [7, [7, p]] + O [A 3 / 2 ] , 

where we have used the scaled phase operator 7 = v^Ay, and we have approximated the 
exponential by its Taylor series expanded in terms of the small parameter yA. After also 
considering the case Ibias < 0, we approximate the Lindblad derived decoherence term in 
the master equation as 

P = [-7, P] ~ [7, [7, P)) + O [A 3/2 ] 

so that in the limit that — > 00 and \/A — > then y/AEb/H is a constant E^/h. In this 
limit the master equation is 

i r ~_ 

which is the washboard potential arising from the bias-current. Thus the tilted washboard 
term arises naturally from our master equation, accompanied by an intrinsically irreversible 
part. 



III. THE QUANTRONIUM MEASUREMENT MODEL 

In order to apply our irreversible current source approach to recent experiments we con- 
sider the Quantronium qubit system depicted in Fig.IU The design of this charge-phase qubit 
is similar to that of the Cooper pair box transistor |23|. The device consists of two identical 
low capacitance Josephson junctions with a coupling energy Ej/2 and capacitance Cj/2. 



S 




FIG. 1: Circuit diagram of the Quantronium qubit. 

These junctions are on either side of the isolated superconducting charge 'island' which is in 
a state of paired electron charge 2eN ', where N is the number of Cooper pairs on the island. 
This island is incorporated into a superconducting loop with a larger Josephson junction, 
which by design, has a coupling energy of E J0 ~ 20Ej and a large shunt capacitance C; 
which was used in the experiment to reduce phase fluctuations. The design of this device 
requires that the characteristic energies Ej and the charging energy Eq = 2e 2 /(Cj + C), 
where 1/C' = 1/C g + 1/4C, are comparable so that neither charge or quantised flux states 
in the loop are good quantum numbers. The discrete energy states of the device are quan- 

nn 

turn superpositions of several charge states[3iJ, |2J|. Control of the qubit is made via the 
pulsed microwave voltage source U (t) which is capacitively coupled to the Cooper pair box 
by the capacitor C g , and the applied flux $ x through the three junction superconducting 
loop. These provide the elementary single qubit manipulations. 

For this device the relation 5 = 7+2e$ x / H between the combined phase 6 = <\>\ — 02 across 
the Josephson junctions of the Cooper pair box, and the phase 7 across the larger Josephson 
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junction provides the readout process of the Quantronium quantum state. From this relation 
the two lowest energy states of the Quantronium have different persistent currents in the 
three junction loop. This difference is used for state readout, a current pulse Iuasif) from the 
'ideal' current source is applied where the height of the pulse is chosen so that the transition 
to a voltage state is made for only one of the Quantronium energy eigenstates; when the 
addition of the loop persistent current state and the current pulse exceeds the critical current 
of the large junction. This process discriminates between the two qubit states associated 
with the two lowest levels of the Quantronium. 

The Hamiltonian for the Quantronium, which we consider in terms of the master equation 
Eq. and its approximation Eq. (J7J), that is without the energy term corresponding to 
the readout current source Iu as is 



Here we have used the terms: the phase operator (p = (</>i + <p2)/2 which is conjugate to the 
Cooper pair number operator N, the dimensionless gate charge N g = C g U/2e, the phase 
bias <p = 27r$ x /$o where $o = h/2e is the flux quantum, and the charge Q on the large 
Josephson junction with shunt capacitance C . In this Hamiltonian we have neglected the 
energy term corresponding to the loop inductance of the device based on the size of the 
device. 

To analyse the measurement induced decoherence in our model we simplify the Hamil- 
tonian H by considering the dynamics of the lowest two qubit eigenstates where we use |0) 
and |1) to denote the lowest and first excited state of the Quantronium system respectively. 
Here we work at the point where the applied flux § x is set to and the dimensionless gate 
charge N g is set to 1/2. In this configuration the qubit energy levels are separated by the 
Josephson Junction coupling energy so we write our Hamiltonian as 



where Eqo = 2e 2 /C and N is the charge operator for the large Josephson junction which is 
conjugate to the phase operator 7. This Hamiltonian describes a two level system separated 




+ 7^ + £jo(l-COS 7 ). 




+ E J0 (1 - cos 7) , 



(9) 
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by an energy Ej(l — cos(7/2)) where the phase 7 provides the coupling between the qubit 
and the readout junction. If we now assume that the large junction is in a localised semi- 
classical state near (7) = 0, then by expanding cos 7 in Eq. Q to second order in 7, we 
obtain the Hamiltonian 

H R = E C0 N> + + (io) 

which has a form of a displaced simple harmonic oscillator. 



IV. THE DC-SQUID MEASUREMENT MODEL 




FIG. 2: Circuit diagram of a DC-SQUID used for qubit state detection via a measurement of the 
qubit's magnetic flux. 

In addition to the Quantronium experiments our approach is applicable to the systems 
where a two level quantum device has been measured by a small inductance DC-SQUID such 
as the persistent current qubit [17]. In these experiments the coherent oscillations in a low 
inductance three Josephson junction qubit structure have been observed. Similarly, the use 

n 

of low inductance microSQUID 32J structures have been proposed to readout the quantum 
state of nanometre scale magnetic particles of large spin and high anisotropy molecular 
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clusters [zJ]. Here the measurement of a magnetic flux quantum state inductively coupled 
to a DC-SQUID with a low inductance relies on the induced change of the effective critical 
current of the the DC-SQUID, for this type of measurement a current ramp scheme is used 
which is similar to that used in the Quantronium readout process. 

In Fig.|2]we consider two Josephson junctions with a coupling strength E J0 /2, capacitance 
Cjo/2 and phases (as shown) of <pi and 02 in a superconducting loop. For this device we 
define the total phase 7 = (0i + <f>2)/2 across the device and the applied flux 

where fia z is the magnetic flux of the qubit state. When the loop inductance is small then 
the flux through the loop $ ~ & x — fi a z , also when the charging energy of the Josephson 
junctions is small so that the quantum state of the DC-SQUID detector is well defined in 
phase and the energy of the first excited state of the detector is larger than the other energies 
of the system so that it exhibits ground state behaviour, then we can write the Hamiltonian 
of the system for the master equation Eq. (jSJ) and its approximation Eq. ((7j) as 



H = H + H : 



s, 

here the DC-SQUID Hamiltonian is 

(eN) 2 

H s = -=7; 1- E Jo (! - cos 7 cos(v?x - 5<p<r z )) 

where tp x — 5(pa z = n(& x — /icr 2 )/$o and the qubit Hamiltonian is Hq = (eo <J Z + to o' x )/2. 
For small 5ip, eliminating the constant terms and assuming that the tunnelling between the 
flux states of the qubit has been turned off, to = 0, we simplify this Hamiltonian H to 

H' = |cr 2 + E co N 2 - E J0 cos (p x cos 7 
- E J0 5ipa z sm(p x cos 7. 

This Hamiltonian is similar to Eq. Q, and since we assume that the DC-SQUID is localised 
near (7) = we again make a second order approximation to the 7 terms to arrive at the 
reduced Hamiltonian 

„ , t2 , Ejo cos (p x 7 2 , Ejo dip sin tp x >y 2 a z 

.Hr = tLCQ ^ H 7, 1 7 5 
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when the qubit energy level separation satisfies eo/2 = Ej 5ip sin^.. Since the form of this 
Hamiltonian is identical to Eq. (fTU|) then the model presented for the Quantronium can be 
directly applied to the measurement of the magnetic flux of a qubit with a low inductance 
DC-SQUID. 



V. THE REVERSIBLE CURRENT SOURCE 

A. The Reversible Current Source Wigner Function 

To investigate the Hamiltonian dynamics of the Quantronium measurement model (and 
by analogy the DC-SQUID measurement model) we consider the simplified Quantronium 
Hamiltonian derived in the previous section: 

H«=E c ,Ni + ^ + ^^. (11) 

We use this Hamiltonian to analyse the measurement induced decoherence relative to the 
washboard potential phenomenology, which does not include the effects of decoherence. In 
this section we derive the decoherence-free dynamics of the system using the standard tilted- 
washboard model by including the term E^j in the Hamiltonian Eq. (fTTj) . In this model, 
the density matrix for the qubit and the readout device evolves according to 

p = -~[H R -E b y,p}. (12) 

We decompose p as 

p= /o + (t) ® |0> <0| + p x (t) <g) [0) <1| 

+ p x f (t)® |l)(0|+p_(t)® |1}(1|, (13) 

where p + and p_ describe the evolution of the Josephson junction when the qubit is in the 
states |0) and |1) whilst p x describes the coherence between them. We assume the initial 
state of the system p(0) is a product state of the readout Josephson junction density matrix 
w and the qubit in the symmetric state (|0) + |l))/\/2, so 

P(0) = ^®(|0> (0| + |0)(1| + |1)(0| + |1)(1|). 

The dynamics of p + and p_ do not depend on p x , so their dynamics are described by the 
qubit-state dependent Hamiltonian 
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We define the two sets of raising and lowering operators a)± and a± by 

7 = r±(c4 + a±) = (a ] ± + a±} 



and 



N 



± 



a± — a± 



(14) 



(15) 



where X±— \j2vj (1 ± p/4), z/ = E C q/E J0 , p = Ej/E J0 , and .ffft± = ftu;±c4_a± where 



2E C0 Ejq I 1 ± 



E, 



AE 



JO 



Using these scalings we can define the two independent equations for p + and p_ 



P+ = -i 
p_ = -i 



uj-a)_a 



E h T + / t 



a+ + a+ j , p 



EbT ~ ( t , , 
— — (ai + a_ J ,p_ 



(16) 
(17) 



By the anti-commutation relation, {A, B} = AB + BA, we define the equation for the 
off-diagonal element p x as 



Px = -* 



a>xa x a x — ,_ 
n 



Eb^x ( t 



< + «x ,Px 



8ft 



o x + o x 1 . 



Px 



(18) 



For the off-diagonal component p x we have defined the raising and lowering operators a x 
and a x where 



7 = T x (a t x +a x ) = YTy (4 + fl x) , 



(19) 
(20) 



A x = v2Ts and fta> x = ^2EcoEj . The master equation for p x defines the dynamics of both 
the off-diagonal elements of the density matrix, where the equation for p x * is the Hermitian 
conjugate of Eq. (JTHJ). 

To solve the dynamics of the system, we transform to a Wigner representation of the 
state [3, 2^1 • To obtain the equation of motion for the Wigner function we first derive the 
characteristic function equation of motion dT(/3)/dt = Tr(Dp), where the characteristic 
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function is defined as T(/3) = Tr(Dp) and D is the displacement operator defined by D = 
exp(/3a^ — 13* a). Writing D in normal and anti-normal order then we can find the relevant 
operator rules for converting to the characteristic function equations. The Wigner function 
equation of motion is found by taking the Fourier transform of the characteristic function 
equation of motion T(j3). Thus 



-oo 

W{a) = I e^ a -^T{f3)d 2 (3. (21) 



After performing this procedure we use a compact notation to write down the Wigner 
function equations from the three master equations Eq. (|16|) - Eq. (II 8J) : for the operators 
a + , a_ and a x defined in the three master equations we correspondingly have the complex 
parameters a+, a_ and a x but we drop the subscripts since they appear separately in the 
three characteristic function equations. This procedure provides us with three uncoupled 
equations 

W + (a) = — iu + (d a *a* — d a a)W + (a) 

- l -^§±{d a -d a *)W + {a), (22) 

W_(a) = - iuj_{d a *a* - d a a)W„(a) 
iE b T 



h 

W x (a) = — iuj x (d a *a* — d a a)W x (a) 
iE b T x 



(d a -d a *)W4a), (23) 



h 



(d a -d a *)W x (a) 



2 



^^(2(a*) 2 + 2a 2 + A\a\ 2 )W x (a) 

oh 

iEjT x 2 



(d 2 a *+cP a -2d a d a *) W x (a). (24) 



We note that each of these three equations are described in three separate co-ordinate 
spaces related to each other by a small scaling factor. This same procedure will be used in 
the description of the irreversible current source described in the following section. Eq. ()24)1 
can be expressed in terms of the phase 7 and charge iV variables using the definitions Eq. 
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(JHJ) and Eq. (JT3J), doing so we find 



W x (7, N) = -u x \X x d 7 N - -^d Nl ) W x (7, N) 

-^-d N W x (^y,N) 
n 



iEj 



(4 7 2 -^)^ x ( 7 ,iV). (25) 



16ft 

B. The Reversible Current Source Wigner Function Solution 

The first two Wigner function equations (Eq. (|2*2*|) and Eq. ([23))) for the readout Joseph- 
son junction density matrix component elements p + and p_ can be solved analytically using 



351 ] since the equations 



the Wang and Uhlenbeck solution for a linear Fokker-Plank equation 
are of the form 

W±(a±, t) = (-V T z .M±.z + V T z .N ± .V z /2) W ± (a ± , t). (26) 

where 

M+ = 




and a = a — EbT±/(fkJ±). From these two solutions W + and W- we can specify the Wigner 
function for the reduced state of the readout junction, since from the definition of the Wigner 
function we have 



/+00 
-00 



=1 (W + + W.) , (27) 

where the trace is performed over the Josephson junction and qubit states. From the Wigner 
function we can obtain a probability distribution for the state of the system in the state 
variables 7 or N by integrating over the state variable for the state variable iV or 7 respec- 
tively. 

This Wigner function for the combined system does not show the coherence that exists 
between the states of the qubit, that is it cannot be used to distinguish between a pure and 
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a mixed state. We therefore construct a function from the three equations for the Wigner 
function terms W + , W- and W x that we derived from the readout Josephson junction 
density matrix component elements p + , p_, and p x ; this function is found by directly Wigner 
transforming both sides of Eq. (|13j) over the Josephson junction degrees of freedom which 
defines the operator 

W.(-y,N,t) = W + I0)<0| + W_|1)<1| 

+ W r x *|l)(0| + W x |0)(l|. 

From this we can calculate the projection onto the initial state W s (7, iV,t) = 
(+\W s {j,N,t)\+) where |+) = + |0)) /y/2 so that 

W s (j,N,t) = l - (W + + W.) + Re(W x ). 

Integrating this function over the canonical coordinates gives the probability to find the 
system in the initial state at time t. 

The solutions for the diagonal Wigner function terms W + and W- obtained from the Eq. 
(|22|) and Eq. (|23|) are the Gaussians 

W ± (a ± ,t) = ^-l^exp ^-l u± T .Cl\u^j , (28) 

where 

Eb£± — iui±t 

1 ^+ 

u±(a ± ,t) 
and the covariance matrix 




e -iu±t \ / e~ iu±t 
C± = I .C . 

e luJ±t I \ e tu±t 

which decays from the initial covariance matrix: 

(a 2 )o-(a) 2 (|a| 2 ) - I («) | 2 

(l«| 2 )o-|(«)o| 2 <(«*) 2 > -(a*)o^ 

The solution (Eq. (|25jl) for the terms W + and W- correspond to Gaussian functions in the 

(a±, a?j_) co-ordinate space. The initial state of the Josephson junction at t = is a Gaussian 

centred about zero, so that at the instant the bias current Iu as is applied 

2 

W s (a x , 0) =- ex P (-2«x«x) 

=^expf-^-A x iV 2 ) , (29) 
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this initial condition implies «o = 0, and the covariance matrix 

C± = C 



I / A x 2 — A± 2 X± 2 + A x 2 



4A±A X I A± 2 + A x 2 A 2 -A± 2 



Lim 

Ejo» b J 




(30) 



From the analytic solutions for the diagonal Wigner function terms W + and W- we see that 
they correspond to fixed width Gaussian curves that rotate on elliptical orbits through the 
(7, N) co-ordinate space at different frequencies uj + and uj- and the centre of the orbits are 
located at E^Y '+/ '(hu + ) and EbT_/(hw_) along the phase axis respectively. These diagonal 
Wigner function terms in the absence of decoherence maintain their width and hence their 
noise characteristics during their evolution. Hence the Wigner function follows a complicated 
periodic motion. For instance, after a certain number of oscillations at time 2T$ = 2n /(uj + — 
ujS) the Wigner function term W + with the larger frequency uj + has completed an extra 
oscillation about its elliptical orbit compared to the Wigner function term W-. 

For the off-diagonal Wigner function term jy x (7,iV, t) we assume a solution of the form 

W x (7, N, t) = exp(a(t) 7 + b(t)N + c(t) 7 2 

+ d(t)jN + e(t)N 2 + f(tj). (31) 

From Eq. (|25|). the coefficients in the exponent evolve according to 

a(t) = ^-b(t) + Gd(t) + (-7 +lE \ {t)d{t), (32) 
A x 2 

b(t) = -X x ua{t) + 2Ge{t) + (-1 + %E) e{t)b{t), (33) 

c(t) = ^d(t) -%E+ ^~ I + %E ^ d{t) 2 , (34) 
A x 4 

d(t) = -2\ x uc(t) + ^e(t) + (-/ + iE) e(t)d(t), (35) 
A x 

e(t) = -X x ud(t) + (-1 + iE)e(t) 2 , (36) 

f{t) = Gb(t) + tl±^l ( 2e (t) + b(t) 2 ) , (37) 

where G = —E^/h, E = Ej/4h and 1 = 0. Using the initial conditions a(0) = 0, 6(0) = 
0, c(0) = -1/A X , d(0) = 0, e(0) = -A x , /(0) = ln(l/vr) we can solve for W x {j,N,t) 
numerically. 
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The plot of the full Wigner function W s (j,N,t) for the readout Josephson junction is 
plotted using the experimental parameters of the Quantronium experiment j^] in Fig. |3J 
The interference fringes, centrally located between the two Gaussians W + and W-, is due 
to the coherence between the qubit states, arising from W x . The number of interference 
fringes present at a particular time t is related to the separation of the two Gaussians in the 
(7, N) co-ordinate space, which increases the further the terms W + and W- are apart. As 
the Gaussians separate the centre of the off-diagonal Wigner term W x follows the trajectory 
shown in Fig. |3J In these figures we see that the main feature of these plots is that the 
noise properties and the interference fringes are conserved over time. In the absence of 
decoherence they continuously evolve with a complicated periodic motion. This will be 
contrasted against the evolution of the state in the presence of the irreversible bias current 
decoherence in the next section. 



w(y,N) 




FIG. 3: Plot of the full Wigner Function W s (j,N,t) for the readout Josephson 
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junction. The 
and they are: 



parameters used in this plot are taken from the Quantronium experiment 
Ej = 0.86k b K, E c = 0.68k b K, E co = 0.0037k b K, E J0 = l8Ak b K and E b = 0.97E J0 The time 
after which the bias current of 0.77 fiA is applied for this plot is t = 16.932 ns. The plot demonstrates 
the different terms that appear in the decoherence free Wigner function for the readout Josephson 
junction. In the plot we can see that at this particular time the two Gaussian curves corresponding 
to the diagonal Wigner function terms W+ and W— are separated and the interference fringes that 
correspond to the off-diagonal coherence term Re(VF x ) appears between them. 
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FIG. 4: Contour plots of the decoherence free Wigner function W s (7, N, t) with the trajectories of 
the centre of the three Wigner function components W+(7, N, t), W-(-y, N, t) and Re(Ty x (7, N, t)) 
superimposed in red, green and blue respectively, (a) The contour plot of the three Wigner function 
terms is for time t = To/4, (b) The contour plot is for time t = Tq. Where the time To is defined 
as the time when the diagonal terms W^.(j,N,t) and W—(j,N,t) are the most separated and sit 
on opposite sides of their respective trajectory ellipse, To = vr/(w + — uiJ). The parameters used in 
this plot are based on those from the Quantronium experiment and they are: Ej = 25 x 0.86k b K, 
E c = O.G8k b K, E C o = 0.0037k b K, T J0 = 18Ak b K and E b = 0.97T J0 ; here the qubit energy Ej has 
been increased by a factor of 25 in order to exaggerate the difference between the trajectories of 
the W+ (7, N, t) and W—(y, N, t) terms. In these plots with the absence of decoherence the Wigner 
functions continuously evolve, moving about their respective trajectories maintaining their height 
and shape and hence conserving the equal probabilities of finding the system in either of the two 
qubit states which is consistent with the qubit symmetric superposition state (|0) + |1})/v2- The 
superimposed trajectories are shown from time t = until t = 2Tq 
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FIG. 5: Parametric plot of the trajectory of the centre of the off-diagonal Wigner function term 
W x (7, N, t) for the parameters used in Fig.^J Here the trajectory is shown from t = to t = 20 x To 
where To = 7r/(u;+ — cj_). From this plot we see the continual coherence of the Wigner function 
under the absence of decoherence, since the off-diagonal W x (7, N, t) term continuously evolves in 
the co-ordinate space without decay. 

VI. THE IRREVERSIBLE CURRENT SOURCE - MEASUREMENT INDUCED 
DECOHERENCE 

A. The Irreversible Current Source Wigner Function Equation 

In the previous section the behaviour of the Quantronium system under the application 
of a bias current, introduced as a Hamiltonian term, was examined. Now we investigate 
the effect of an irreversible bias current model by adding the Lindblad derived terms that 
appear in Eq. (JJJ) to show the relative decoherence in the system's evolution. In this case 
we have the system density matrix p defined by the master equation 



Using the scaling factors Eq. (|14|) and Eq. (|19|) we can write the following three master 
equations to describe the elements of the density matrix p as 



[H R -E bl ,p] 



E b 



[7, h,p}}- 



p 



2h 




E b \T 

2h 



+ 



a' + + a + , a' + + a + ,p + 
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P- 



— I 



t E t> T - 

u;_a_a_ — 

n 



at + a_ I , p 



\E h \TJ 
2h 



at + a_ 



a_ + a_, p_ 



— i 



u; x a x a x 



E b T x f | 



I -^6 1 r x 
2h 

iEjTy 2 
8h 



a x + a, x , 



at + a 



Ox +°x J ,P 



a f x + a x ,p x 



x ' x / , yx 



Following the same procedure used in Section we obtain the three component Wigner 
function equations. The three uncoupled Wigner function term equations are 



W + (a) = — icu + (d a *a* — d a a)W + (a) 



Wja) 



\E b \T + 

2h 

iE b T + 

h 



-{dl+d 2 a -2d a d a *)W+ 
(d a -d a *)W + (a), 



.a 



iuj-(d a *a* — d a a)W-(a) 
\E h \T_ 2 



2h 
iE b T_ 

h 



-{d 2 a *+dl-2d a d a *)W_i 
(d a -d a *)W4a), 



a 



(38) 



(39) 



W x (j, N)=-lu x {\ x d y N - —dNlj W x (r/, N) 

- ^d N W x ( 7 , N) + l -^d 2 N W x ( 7 , N) 
iEj 



16ft 



(4 7 2 -^)W x (7,A0, 



(40) 



where the simplified a notation convention from the previous section has again been used. 
By solving these three equations we can investigate the evolution of the Quantronium device 
in the presence of the irreversible bias current and in particular the decay of the off-diagonal 
term W x , that projects the qubit into one of its eigenstates with probabilities related to the 
initial state of the qubit. 
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B. The Irreversible Current Source Wigner Function Solution 



As was the case for the Hamiltonian evolution of the Quantronium system with the 
reversible bias current term, which we described in the previous section, the first two Wigner 
function equations (Eq. (|55j) and Eq. (J2HJ)) can be solved using the Wang and Uhlenbeck 
solution for a linear Fokker-Plank equation since the equations are in the form of Eq. ()26|) 
where 



M. 



± 




N- 



± 





' -\E b \r ± 2 


+\E b \r± 2 \ 




V +l^|r± 2 


-\E b \r ± 2 ) 



da\ _ (OL 

da*)' \a* 

and a = a — E b T±/ {ftu)±). The solution for the diagonal Wigner function terms W + and 
W_ is the again the Gaussian 

W ± {a ± ,t) = ^-l^exp ^ u± T .Cl\u^j , (41) 

where 

E b T± _ „-iw±t ( ' _ E b T± 



and from the initial condition Eq. (J2U|) we have «o = and the covariance matrix 



|g fc |r± 2 t »|g 6 |r ± 2 / 2to;±t _ 1 N 

2H 4Hlu± ^ e 



+ .Co- (42) 

e +lu±t ) \ e +luJ±t 

which decays from the initial covariance matrix Co given by Eq. ()30j) . The solutions W + 
and W_ (and W x ) are shown in Fig. |Ul the parameters used in this figure are the same as 
those used in Fig. 0] and Fig. El where the parameters of the Quantronium experiment have 
been used to demonstrate the evolution of the Wigner function with the exception that the 
qubit energy Ej has been increased by a factor of 25 to exaggerate the separation of the 
states in the presence of the increasing noise characteristics of the irreversible bias current. 
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In Fig. IHlwe see that the trajectories of the diagonal Wigner terms W+ and W- through 
the (7, N) co-ordinate space are identical to those found using the reversible bias current 
approach of the previous section. With the application of the bias current Iu as the two 
Gaussians start from the initial condition where they superimposed on each other at the 
origin and then W+ and W_ separate as they begin to rotate about the points E b T + /Huj + 
and EbT-fhuj- on the phase axis with a frequency uo + and uo^ respectively. However, the 
shape and hence the noise characteristics of the W+ and W- terms have changed. The 
off-diagonal elements \E b \T± 2 t/2h in the covariance matrix Eq. (|42jl mean that during the 
evolution of the states, energy from the system is 'leaking' and causing the Gaussians to 
become broader as they separate. For long time this means that the states become virtually 
indistinguishable in the (7, N) co-ordinate space. 

For the off-diagonal term W x (7, N) we solve Eq. (|40|) with a solution in the non-positive 
definite form Eq. (|31|). Using this form of solution we can derive the set of six coupled 
differential equations Eq. (|52jl - Eq. pTJ) where G = -E b /h, I = -2\E b \/H and E = Ej/Ah. 
We can solve this set of equations numerically using the same initial conditions used in the 
previous section, i.e. Eq. (J2HJ). In Fig. El we see that W x decays as the system evolves, 
so that by the time the diagonal terms W + and W_ are the most separated at time T Q it 
has virtually decayed to zero relative to the diffusing and larger diagonal Wigner terms. 
From this numerical solution of the off-diagonal Wigner function term W x we are now in 
a position where we can examine the time it takes for the coherence of the initial qubit 
symmetric superposition state to be lost. Also the effect that this description of the bias 
current as a Poisson distributed kick process has on the qubit when white noise in current 
source is considered. 

VII. DECOHERENCE IN THE QUANTRONIUM EXPERIMENT 
A. Dephasing Time of Read-Out Current Pulse 

In Section|Sand Section IVll we have numerically obtained the solution for the off-diagonal 
Wigner function term W x both with and without the presence of decoherence from the 
irreversible bias current. These numerical solutions were obtained in terms of the functions 
a(t), b(t), c(t), d(t), e(t) and f(t) that specify the off-diagonal wigner function in the form 
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FIG. 6: Contour plots of the Wigner function W s (^f, N,t) in the presence of the decoherence 
from the irreversible bias current, (a) Contour plot at time t = To/4, (b) Contour plot at time 
t = Tq where To = tt/(u;+ — ujJ). The superimposed trajectories show the time evolution of the 
centre of the three Wigner function components W+("/, N, t), W_(7, N,t) and Re(W x (7, N, t)) 
shown in red, green and blue respectively. The parameters used in this plot are identical to 
Fig. 0] and derived from the Quantronium experiment and again the qubit energy Ej has been 
increased by a factor of 25 in order to exaggerate the difference between the trajectories of the 
W+(7, N, t) and W-(~f, N, t) terms. In these plots we note that the role of the irreversible bias 
current decoherence on the Wigner functions. The qubit state which is initially the symmetric 
superposition state (|0) + |l))/v / 2 evolves to a state where the coherence term Wxij, N, t) decays 
to zero and corresponds to a classical equal probability mixture of the two states |0) and |1). 
During the decay of the off-diagonal Wigner function term W x (7, N, t) the trajectory of its centre 
spirals towards a central point as the coherence of the qubit state is lost. As the coherence is 
lost the diagonal Wigner function terms W+^j, N, t) and W-i^y, N, t)broaden out and become less 
localised. The superimposed trajectories are shown from time t = until t = 2Tq. 
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of Eq. (|31|). From these solutions we can determine the role of the irreversible bias current 
on the coherence time of the qubit when it is initially prepared in a symmetric superposition 
state. We calculate the coherence time of the qubit from the length of the Bloch vector B(t) 
for the state of the system defined by 



B(t) = \/{cy x ) 2 + {a y ) 2 + {a z y 



For our coupled readout Josephson Junction and Qubit system we can write the expectation 
of some operator A which operates on the Qubit as 

(A)=Tr(Ap) = Tr n (Tr Qu (Ap)). 



Since the Wigner quasi-probability distribution function 33] allows us to compute expecta- 
tions of operators straightforwardly, that is, 

/+oo p+oo 
/ A(7,JV)W(7, N)djdN 
00 J —00 

where A (7, N) is the Wigner transform of the operator A, we therefore have 

/ d-ydN(W+{0\A\0) + W-(1\A\1) 
00 J —00 

+ W X (0\A\1) + W X *(1\A\0)) (43) 

for some operator A acting on the qubit. From this expression we can obtain the expectation 
values of the Pauli matrices and they are 

/+00 r+00 
/ Re(W x )d^dN 
-00 J —00 

(<t v )= / / 1m(W x )d>ydN, 



00 j —00 



and 

1 p+00 r+oc 



1 /*+oo r-\-oo 

- / (W + - W_) d-fdN = 0. 

2 J — 00 J — 00 



2 

From these expectation values the length of the Bloch Vector for the qubit can be written 

as 



B(t) 



00 r+oD 



—00 j —00 



(44) 



W x (-f,N,t)d-fdN . 

Since W / x (7, ^V") is in the form Eq. (|31|) we can integrate this analytically and then find B(t) 
using the numerical results for the functions a(t), b(t), c(t), d(t), e(t) and fit), doing so we 
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exp [/(*)] 



have 

Bit) = . 2?r 

V4c(t)e(t) - d(t) 2 

\b{t) 2 c{t) - a{t)b{t)d{t) + a{t) 2 e(t) 
X 6XP [ d(t) 2 - Ac(t)e(t) 

In Fig. [7| we show the qubit Bloch vector length as it evolves in time for both the de- 
coherence free and irreversible bias current solutions. In these plots the parameters of the 
Quantronium experiment have been used, and the graph shows the system evolution from an 
initial qubit symmetric superposition state when the manipulation of the qubit has ceased 
and the time scale starts at the instant the read-out bias current pulse is applied. 

The main feature of the decoherence free plots is the periodic nature of B(t), in the 
absence of decoherence the state of the qubit evolves from an initial pure state to a mixed 
state and then back to a pure state when the two diagonal Wigner function terms W + and 
W- are superimposed on each other at time 2Tq. With the application of the irreversible 
bias current to the Quantronium experiment we see that the state's progression to a mixed 
state is hastened and there is no revival of the qubit state back to a pure state at time 
2Tq. From this plot we can see that the Bloch vector length is 0.5 at time 0.034To, or 0.18 
nanoseconds, after the bias current read-out pulse has been applied to the system. 

In the Quantronium experiment the readout pulse lasts for a duration of the order of 0.1 
microseconds, meaning that according to our irreversible bias current model the state has 
been dephased on a time scale that is around a thousand times faster, before any classical 
information about the state has been returned to the experimentalist. The consequences of 
this is that in the experiment the qubit decoheres much faster than the time taken for the 
measurement. 



B. Decoherence Due to Thermal Fluctuations in the Current Source 

The irreversible bias current model presented so far not only has a dephasing effect when 
the readout current pulse is applied to the Quantronium circuit but also during the presence 
of noise in the current source. Here we look at the decoherence that our model predicts 
during the qubit manipulation stage when the mean of the readout current source is zero 
but has white noise fluctuations due to thermal noise in the resistor network which is used 
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FIG. 7: Plots of the qubit state Bloch vector length B(t) versus time for the decoherence free case 
(shown in red) and according to the irreversible bias current model (shown in blue). The time 
scale is parameterised in terms of the diagonal Wigner function term maximum separation time 
Tq. Here we note that in the case of the decoherence free plot the qubit state is a pure state at 
time 2uTq where n is a positive integer. 



in conjunction with t 
readout current pulse 



le voltage source in the Quantronium experiment to implement the 



3l|. In this regime we consider the full master equation that includes 



the irreversible bias current decoherence term 



P = 4^^ + ! f [7,P]-^[7,[7,P]] 



and we write this in the operator form 

p = C H {p} + e a C a {p} + t b C b {p} (45) 

where 

Ch{p} = ~[H iP ], (46) 
eaC a { P } = tf[j,p], (47) 



and 

fi.C.1. { n\ = — 

2h 
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e b C b {p} = —§-[7, [l,p]\- 



Here e a and e& are small parameters compared to the qubit Hamiltonian H and they are 
determined by the zero-bias current noise where e a = Eb and e b = \Eb\. We assume Ef, 
fluctuates due to thermal noise in the external circuit, having zero mean. Therefore Ej, 
denotes the time fluctuating and zero mean current noise. For statistical purposes we treat 
e a and eb as independent variables. Since (e&) = (\E b \) ^ 0, we ignore fluctuations in \E b \, 
whilst (e a ) = (Eb) = and so we consider the effect of thermal fluctuations in Eb- We write 
the solution to Eq. Po"]) as a correction to the exact solution po of the noise free master 
equation 

i 

po = ~T [H,po\ , 

that is, we can write the solution in the form p = po + e a p a + ebPb + 0(e 2 ). That is, p a and pb 
represent the change of the density matrix due to the effect of fluctuations in the coherent 
term and for the intrinsic dissipative term (that depends on \E b \) respectively. Substituting 
this solution p into Eq. (jlo]) and expanding the master equation to 0(e a ) and 0(e&) gives 

Po = £hPo, (48) 

p a = C H pa + £aPo, (49) 



pb = C H p b + C b po- (50) 
Defining p a = p + t a p a and p b = Po + t b p b we find 

Pa = ~[H,p a ] + l -^[l,p a ] (51) 

and 

h^~[H,p b ]- l -^h,h,P b }} (52) 

This allows us to treat the intrinsic dephasing of the current source and the dephasing due 
to thermal fluctuations in the washboard potential independently. Since Eb fluctuates, the 
second term of Eq. (J51|) results in extra decoherence, on top of the intrinsic decoherence 
due to current passing through the readout Josephson junction described by Eq. (p)2*|) . From 
these equations we can estimate the magnitude arising from these different effects. 

Assuming the thermal fluctuations of Eb are well approximated by white noise with zero 
mean we derive a master equation from Eq. (}5T|) which describes the dephasing effect of 
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a fluctuating current. By integrating Eq. ()51|) and substituting into the original master 
equation we have 



Pa = ~ ^ [H, Pa] 



+ 



h 



i r 

E b (t)^,--J ds[H - E b (s)^,p a (s) 



(53) 



Taking the ensemble average of this equation and using (Eb(t)) = and (E b (t)E b (s)p(s)) 
(E b (t)E b (s))(p(s)) (since (pit)) is independent of future noise fluctuations) then 

i 



(Pa) 



[H, (Pa)} 



h 

1 « 



The noise correlation function satisfies 

(E b (t)E b (s)} 



ds (E b (t)E b (s)) (p a (s) 



2tt 



(hias(t)I b ia S {s)) 



5, 2 (0)5(t- S ), 



where by definition 



1 

2^ 



e wt (I bias (t)I bias (0))dt 



is the noise spectrum of the current fluctuations due to thermal noise. Thus, in the presence 
of noise, the master equation is of the form of Eq. (|52j) . That is, it can be written as 



(Pa) = ~ [H, (Pa)} 



2nh 



[7, [7, (Pa)}} ■ 



Now for the Quantronium circuit with a current source output resistance Ri at temperature 
T and an effective input resistance of R2, the current noise spectrum in terms of the thermal 

371 ] Sy(uj) is 



voltage noise spectrum 



Sj(0) 



S v (0) y/AR^hT 



Ri 



(54) 



The ensemble average of Eq. (J52)) is 



(p b ) = ~ [h, (m - [7 ' {m ' 

which establishes that the two kinds of dephasing have the same form. The rates due to 
the intrinsic dephasing of the current source and noise in the washboard potential are given 
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by ^deph = (\Ibias\)®o/(4:irh) and ^ noise = (SV(0)$ /(27rft)) 2 respectively. We have so far 
assumed that Ef, is 5-correlated. Under this white noise assumption (\Ef,\) is singular so 
instead we estimate it from /rms- The thermal noise in the current source is a result of the 
random scattering of electrons in the output resistance and this produces a statistical dis- 
tribution of phonons in the phonon modes of the resistor. In the thermal state each phonon 
mode has a Gaussian distribution for the resultant voltage and hence current fluctuations. 
Summing over the phonon mode distributions in the resistor we obtain the total current 
fluctuations which is also Gaussian distributed. For Gaussian processes 

[2 

(\hias\) = A / —IRMS- 



and therefore we have 



y/ARxfoTB 

'RMS 



\J (hias 2 ) 



R2 

where B is the bandwidth of the Quantronium circuit. 

We can now compare the sizes of the two dephasing terms by referring to the details of 



the Quantronium experiment 



3 it ] . Analysing the readout circuit we can see that the thermal 



noise is produced by a 10kf2 resistor in series and a 50f2 resistor in parallel with an ideal 
voltage source at the temperature of the helium bath. Also the thermal noise from these 
two resistors contribute to the fluctuating current that flows into the Quantronium Circuit 
via a 3.5 kQ input resistance which has a 200 MHz bandwidth. From these parameters we 
are able to determine that 

7™=(^^Y = 9.68 GHz 



and 



(\hias\)®0 ^RMS^O ,rr MU 

Ideph = -j-z = —7=r = 555 MHz, 

471-/1 V87T% 



meaning that the dephasing rate intrinsic to the irreversible bias current is about 20 times 
slower than the rate due to fluctuations in the titled washboard potential. From the relative 
scale of these two terms we can see that the dephasing during the qubit operation will 
be dominated by the fluctuations in the washboard potential, rather than the intrinsic 
irreversible bias current induced dephasing. However, we note that the introduction of the 
irreversible current source into the modelling process has still provided a dephasing effect 
of considerable size relative to the effect of thermal noise in the current source in this 
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case, and therefore may be important for the consideration of other similar current biased 
superconducting circuit experimental models. 

VIII. DISCUSSION AND CONCLUSION 

In this paper we have analysed the bias current readout process of the superconducting 
qubit structures such as the Quantronium (and by analogy those qubits whose quantum state 
is measured by a DC-SQUID like the persistent current qubit). By introducing an irreversible 
bias current term through Lindblad operators that describe the addition and subtraction of 
electrons across the readout Josephson junction, at a rate given by the bias current, we are 
able to obtain a master equation that can be approximated to first order by the Hamiltonian 
washboard potential model - a model that is used throughout the superconducting quantum 
device literature. Therefore this master equation incorporates an additional term to the 
washboard potential terms that dictate the decoherence of the qubit through its coupling to 
the readout Josephson junction. 

The decoherence is a result of the bias current 'counting' the number of electrons that 
pass through the measurement Josephson junction. We propose that such an effect is in- 
trinsic to the application of the bias current to the system and has a cumulative effect of 
decohering the system as electrons pass through the readout Josephson junction. By approx- 
imating the Hamiltonian terms by a harmonic oscillator coupled to a qubit in the symmetric 
superposition state were are able to analyse the measurement induced decoherence before a 
tunnelling process out of the washboard potential occurs and produces a measurable voltage 
for the experimentalist. Looking at this model in terms of the Quantronium experiment we 
have been able to construct the Wigner function for the Josephson junction and analyse the 
dephasing effect upon the application of the external bias current. 

By analysing the Quantronium system we have found that the effect of describing the 
readout bias current in terms of the Lindblad operators is to produce a qubit dephasing time 
of 0.2 nanoseconds after the bias current has been applied. In the Quantronium experiment 
the bias current pulse was applied for a duration of the order of 0.1 microseconds, meaning 
that the state of the qubit has been reduced to a mixed state before the the tunnelling event 
from the washboard potential is observed. Our model changes the understanding of the 
measurement process of the Quantronium qubit and it means that the point of measurement 
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is not the tunnelling event out of the washboard potential but instead arises as a consequence 
of coupling a current biased Josephson junction to the qubit state. Additionally this model 
does not produce extra sensitivity to noise in the current source since by adding thermal 
noise to the irreversible bias current model we showed that thermal noise in the washboard 
potential produces the dominant dephasing effect by an order of magnitude. 

Experimental validation of our model could be predicted by using small current pulses 
during the Ramsey fringe experiment demonstrated by Vion et al. ^| since the role of 
the irreversible bias current is to dephase the qubit. Small amplitude and short duration 
current pulses could be applied to the Quantronium between 7r/2 pulses of a Ramsey fringe 
experiment. Our model would predict that for larger current and longer duration pulses 
the dephasing would become larger and hence influence the decay time seen in the Ramsey 
fringes. In addition to this the process shown in this paper of adding the irreversible bias 
current to the current biased Josephson junction qubitsjia. \i\ could be employed to look 
at the effect of the constant current through the Josephson junction and its resulting deco- 
herence. In this case due to the utilisation of excited states of the washboard potential for 
the qubit states and readout, an appropriate replacement to the harmonic oscillator used in 
this paper would need to be employed. 
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APPENDIX A: ANALYTICAL CALCULATION OF THE OFF-DIAGONAL 
WIGNER FUNCTION WITH REVERSIBLE CURRENT SOURCE 

Using the commutation relation [j,N] —i which gives e~ ldN \ , y) = Ij+O) and the definition 
of the Wigner function 
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we can write 

9 r+oo 

W^, N) = - e^ N (7 - W2|p|7 + rh/2) d Vx . (Al) 

J-00 

Using this form we can calculate the off-diagonal Wigner function VT X (7, iV) that corresponds 
to the density matrix component p x (t) when we decompose the combined qubit and detector 
density matrix into the form: 

P (t) = P+ {t)\o) (o| + P _(0|i)(i| 

+ p x (t)|0)(l|+p x t (t)|l) (0|. 

We are able to calculate the off-diagonal Wigner function W x (7, N) in terms of the wave- 
functions IV'o(O) an< ^ IV'i(^)) which correspond to the single mode Gaussian wavefunctions 
of the Hamiltonian Hr in the qubit eigenstates |0) and |1) respectively since 

A>x(*) = llM*)XlM*)l- 

The wavefunctions IV'o(O) anc ^ IV'i(O) evolve according to the Hamiltonians H R+ and 
H R _ respectively where 

Hr± = fvjja)a ± hx (a + c^) 2 + fre (a + aj . 

The 7 space wave function for the most general single mode Gaussian pure state is 

( 7 |VW> = (27r<(A 7 ) 2 » _1 exp^ 7 /2)exp(^iVo7o/2) 

x exp(iiV 7) exp(-a (7 - 7o ) 2 /2) (A2) 

where 

7o = (7) , 
No = (N) , 
a — <7i + i<7 2 , 

and 

(A 7 AAQ sym (jV 7 ) + (7iV)-2(7)(iV) 
2<(A 7 ) 2 > 2((A 7 ) 2 ) 
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The phase angle <5 7 is set to zero. By using the single mode Gaussian form Eq. (jA2|) and 
the Wigner function definition Eq. (jAl|) we can calculate the off-diagonal Wigner function 

2 r+ao 



2 r+oo 



rr .l-oo \ 2 



- J-00 \ 2 



4>o (t))(Mt) 



7 + 7f)*fe 

7 + y } dr) x . 



This integral is in the form 



where 



and 



C 



1 e 2 

- GXP 2A 



B = iN--{(N) + + (N)_) 
+ ^ ff +(7 - (7)+) - 2 - cr -(7 - (7)-), 



A = ~o\ 
4 4 



1 * 
-CJ* 

4 " 



C 



^<(A 7 ) 2 > + ((A 7 ) 2 >. 



7T 



x exp I ~«iV) + < 7 ) + - (iV)_< 7 )_) + z(iV) +7 - *(iV>_ 7 

x exp f-^+(7 - (7>+) 2 - - (7>-) 2 

Here we have used the notation ( 7 )±, (N)±, ((A 7 ) 2 ) ± and a ± to distinguish the mean and 
noise parameters of the single mode Gaussian states \ipo{t)) and \ipi(t)) respectively. 

In order to fully specify the off-diagonal Wigner function we need to calculate the quan- 
tities (7), (N), a and S x for both the states \ipo(t)) and \ipi(t)) where 



(7) 



«a> + <« T » 



(TV) 



A 



((«) " <a f » 



((A 7 ) 2 ) =-{{aa) + (aV) + (aa f ) + (a f a) 
-2(a)(a t ) - (a) 2 - (a f ) 2 ) 
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and 



(A 7 AiV) sym =(iV 7 ) + ( 7 jV> - 2( 7 >(A0 

= l -((aW)-{aa) + (a)(a)-(aW))- 

To calculate these quantities we find the two sets of equations that solve for (aa), (a^a^), 
(at a), (aa 1 ), (a) and (a)) for the two qubit-eigenstate Hamiltonians for the system qubit and 
detector H R ± in the qubit states |0) and |1) respectively. In order to find these we construct 
the set of six Heisenberg equations of motion for each Hamiltonian using the relation 

dA/dt = -i[A,H R ±]/h (A3) 

and solve them simultaneously. From our Hamiltonians H R+ and H R _ respectively we find 
the set of two coupled differential equations 

d(a)/dt = - iu(a) =F 2ix{(a*) + (a)) - it 
d(a J< )/dt = + iu(a ] ) ± 2ix((a j ) + (a)) + ie, 

which we solve using the initial conditions (a(0)) = and (a^(0)) = 0. The four remaining, 
coupled equations of motion for H R+ and H R _ are 

d{aa)/dt = — 2i (uj ± 2%) (aa) 

=F 2ix{(a)a) + (aa))) — 2ie(a) 
d(a)a))/dt = + 2i(uj± 2 X ) (aW) 

± 2ix((a ] a) + (aa))) + 2ie(a)) 
d(a)a)/dt = =p 2ix((a)a)) — (aa)) — ie({a') — (a)) 
d(aa))/dt = q= 2ix((a)a)) - (aa)) - ie((a ] ) - (a)), 

Using the solutions (a)t and (a')t we write the four coupled differential equations in matrix 
form V t = A.x + v, where x = (aa, a) a* , a) a, aa) , a) , a) T , V t contains the time derivatives 
of the components of x, and v contains the terms containing (a) t and (a)) t . We solve this 
system of equations by diagonalising the matrix A by forming the matrix D containing its 
eigenvectors corresponding to the eigenvalues 

jo, 0, +2iy/u 2 ±4ux, -2i^u 2 ±AuJx} ■ 
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Once in the diagonal form D _1 .V t = D^A.D.D^.x+D^.v we can solve the four uncoupled 
differential equations and then transform back the solution to the original basis. Using the 
initial conditions (aa) = 0, (aW)o = 0, (a^a)o = 0, (aat)o = and (aJa)o = 1 we have the 
solution for 7 , N and o: 



v / 2Ae(cos {t^/oo 2 ± Auox) ~ 1) 
y/2e sin [t^u 2 ± Aux) 



(N) 



± 



a± = 



y/X(cU 2 ±Au X ) 

uj 2 ± 2 X {2uj + i^/uJ 2 ±4uj X sm (2t^uo 2 ± Auox)) 



<(A 7 ) 5 



± 2x(l + cos {2t^/u 2 ±Aux))) 
\{u ± 2x(l + cos (2t^uj 2 ±AuJx))) 



>± 



2(cu±4 X ) 

In this calculation we have set the phase angle 8 1 to zero for both the IV'o(^)) and 
states. Now that we have fully specified the mean and noise parameters for the two single 
mode Gaussian sates \ipo{t)) an d IV'i(O) we can write this solution in the form W x { r f ) N) = 
exp(a(t)7 + b(t)N + c(t)7 2 + d{t)^N + e(t)N 2 + f(t)) where 

(i(N) + + i(N)_ + a + ( 7 >+ ~ ~ * + ) 



a(t) 



+i(N) + -i(N)_ + a + { 1 ) + + a*_( 1 )_ 



b(t) 



2((JV) + + (JV)_)-2t - + (7) + + 2» ff !( 7 ). 
<j + + <j*_ 



r* "\2 



w 2 + 2 - 2(a + + a*) 



d(t) 



2ia* — 2ia 4 



e(t) = 



a, + a* 
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and 

/(*) = ln I , 

+\{(N) + (l) + -{N)-(l)-) 

-\°+(l)+ 2 ~ \ a *-^)- 2 
(i(N) + + i(N)_ + <r + (>y) + -<T*(>y)_) 2 
2(a + + a*) 

APPENDIX B: ANALYTICAL CALCULATION OF THE OFF-DIAGONAL 
WIGNER FUNCTION WITH IRREVERSIBLE CURRENT SOURCE 

For the off-diagonal term W / x (7,A r ) of the Wigner function including measurement in- 
duced decoherence we solve (Eq. (}4*0|) ) with a solution in the non-positive definite form 

W x (7, N, t) = exp(a(t) 7 + b(t)N + c(t) 7 2 

+ d(t) 1 N + e(t)N 2 + f(t)). 

Using this form of solution we can derive the set of six coupled differential equations: 

d(t) = ^-b(t) + Gd(t) + ( ~ J r ~!~ lE \ {t)d{t), 

A x 2 
b(t) = -X x ua{t) + 2Ge(t) + {-I + iE) e(t)b(t), 

c(t) = —d{t) -iE + ± J -d{t) 2 , 

X x 4 

d{t) = -2X x uc(t) + -r^e(t) + {-I + iE) e{t)d(t), 

X x 

e{t) = -X x ud(t) + (-1 + iE) e(t) 2 , 

f{t) = Gb(t) + tL+13 (2e(t) + b{tf) , 

where G = —E^/h, I = —2\Ef,\/h and E = Ej/Ah. This system of equations is solved 
by first considering the three coupled equations for c(t), d{t) and e(t) whereby using the 
transformation of variables 



4i 

c 

Ai 



d~^{i~-Ke) 
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and 

(t + A x g ) 

u= — 4 — 

we have 

z = 2iu x z - (-1 + iE) (z - uf /2 - 2iE, 
z = -2iu x 'z + (-1 + iE) (z + uf /2 + 2iE, 
u=-(-I + iE) (z -u)(z + u) /2 - 2%E. 

Using the relation dZ/dt = -(-I + iE)Z(z - z - 2u), where Z = zz + u 2 - iE/(-I + iE), 
we use a second transformation of variables 

U = u/Z, 
A = z/Z 



and 



so that 



A = z/Z 



JJ = _(-/ + iE) {zz + u 2 + ME/ {-I + iE)) 



A = 2iuj x A 
A = -2iuo x A + 



2 (z'z + u 2 -ME/(-I + iE)Y 
(-1 + iE) (zz + u 2 + ME I (-1 + iE)) 
2 (z'z + u 2 -AiE/(-I + iE)Y 
-I + iE) (zz + u 2 + ME I (-1 + iE)) 



2 (z'z + u 2 -ME/ (-1 + iE)) 
From these equations we construct the differential equation 

W + 4 ^ x S~ " 16c ^- J + lE ^ p = 

where 

p _ (z'z + u 2 + ME I (-1 + iE)) 
~ (z'z + u 2 -ME I (-1 + iE)) 

The solution P is the sum of exponentials 

P = C ie Xlt + C 2 e~ Xlt + C 3 e X2t + C 4 e~ X2t 

where 



Ai,2 = ^J-2cu 2 x ± 2uj x ^Juj 2 x - ME 2 - ME I 
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and 



Ci 



A 2 2 (2J - i£ ± P(P + 4i/)/(2Ai) 



2 ,1 r/ \ 2 \ 2 



4/(A 2 2 - Ai ) 



Ai i! (2/-i£?±£?(f? + 4i/)/(2A 2 ) 
C3 ' 4 " 4/(A 2 2 - Ai 2 ) 

From P we have the solutions for c(t), d(t) and e(t) 

c(t) 

d(t) = 
e(t) 



2 dP | d^P 



4A x w 2 (-/ + iP)(l -P) 

1 rf 2 P 
~2u?(-/ + iP)(l -P)ltP 
A x rfP 



2(-/ + iP)(l -P) dt' 

If / = 0, Ai j2 are pure imaginary and so P simply oscillates resulting in oscillatory solutions 
for c{t), d(t), e(t). The remaining coeffcients a(t) and b(t) are coupled together, satisfying 
a forced, parametrically excited second order ordinary differential equation. To see this let 



m d(y(t)y/l-P(t)) 
a(t) = ; !L - 

2(1-P(t)) v /2A7 dt 

m = JK ( m + 4G 



2^1 -Pit) (-I + iE) 
Then y(t) satifies the following equation. 

#y ( 1 ( 1 dP(t) V\ 4iGy/T=PV) 
dt 2 ^y 2\(1-P(t)) dt J ) yy > -I + iE 

Since initially both a(0) = and 6(0) = their subsequent solution is proportional to G 
and 

b 2 {t)c{t) - a{t)b{t)d{t) +a 2 {t)e{t) 
4c(f)e(f) - d 2 (t) 

is proportional to G 2 . Initially this is zero and for small times is quadratic in time. The 
solution to f(t) is found through integration. 
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